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Abstract 

Recent insights to melt migration beneath ridges suggest that channelized flow is a 
consequence of melting of a heterogeneous mantle, and that spreading rate modulates 
the dynamics of the localized flow. A corollary of this finding is that both mantle het¬ 
erogeneity and spreading rate have implications for the speed and time scale of melt 
migration. Here, we investigate these implications using numerical models of magma 
flow in heterogeneous mantle beneath spreading plates. The models predict that a broad 
distribution of magma flow speeds is characteristic of the sub-ridge mantle. Within the 
melting region, magmatic flow is fastest in regions of average fusibility; surprisingly, 
magmas from sources of above-average fusibility travel to the ridge in a longer time. 
Spreading rate has comparatively simple consequences, mainly resulting in faster seg¬ 
regation speeds at higher spreading rates. The computed time scales are short enough 
to preserve deep origin ^^°Th disequilibria and, under favourable parameter regimes, 
also ^^®Ra excesses. An important prediction from the models is that mantle hetero¬ 
geneity induces significant natural variability into U-series disequilibria, complicating 
the identification of relationships between disequilibria and ridge properties or chemical 
signatures of heterogeneity. 

Keywords: U-series disequilibria, mid-ocean ridge basalt, mantle heterogeneity, melt 


* Corresponding author 

Email addresses: smw@geus.dk (Samuel M. Weatherley ), richard.katz@earth.ox.ac.uk 
(Richard F. Katz) 

^Present address: GEUS, Geological Survey of Denmark and Greenland, 0ster Voldgade 10, 1315 
Copenhagen K. Denmark. Tel: -1-45 91 33 38 68 


Preprint submitted to Geochimica et Cosmochimica Acta 


September 8, 2015 



migration, mantle melting, pyroxenite 


1. Introduction 


Chemical heterogeneity in Earth’s mantle induces local and regional perturbations 
in the fusibility and melting rate beneath oceanic volcanoes. For example, regions en¬ 
riched in subducted oceanic lithosphere start to melt at deeper depths and with higher 
melt productivity than mantle peridotite (Hirschmann and Stolper, [1996 ), whereas 
other regions are so refractory, and contribute so little to the genesis of magmas, that 


they are essentially unsampled by mid-ocean ridge basalts (Liu et ah, 2008). In two 


recent numerical studies Weatherley and Katz (2012) and Katz and Weatherley (2012) 


considered the effect that mantle heterogeneity has on the dynamics of melt migra¬ 
tion beneath oceanic volcanoes. Their models suggest that (i) heterogeneity causes 
magma to localize into a network of high porosity channels, and (ii) that the topology 
of the channel network depends on the spatial arrangement of mantle heterogeneity. 
Channelised flow is often invoked to explain the major and trace element character¬ 


istics of oceanic basalts and mantle rocks (e.g. Kelemen et ah, 1995, 1997), and to 
explain geochemical evidence for rapid magma flow speeds in the mantle (see reviews 


by Lundstrom et ah, 2005 Elliott and Spiegelman, 2014). In this paper we explore the 


prediction arising from Weatherley and Katz (2012) and Katz and Weatherley (2012) 


that mantle heterogeneity has important but as yet unrecognized consequences for the 
time scales of melt migration beneath oceanic volcanoes. 

Temporal constraints on magma flow in the mantle largely are derived from observa¬ 
tions of U-series disequilibria in very young or zero age MORE. Specihcally, disequilib- 
ria provide temporal constraints on combined melt generation and transport processes 
that have occurred within approximately 5 half-lives of the shorted-lived nuclide in a 
parent-daughter nuclide pair. For understanding magma flow processes, the most useful 
nuclides are ^^°Th, ^^^Pa, ^^®Ra and ^^°Pb; these decay with half-lives of c. 75 ka, 33 ka, 
1.6 ka and 22 yr. Since a wide range of processes can generate U-series disequilbria, 
models are often used to explore the implications of disequilibria for the time scale of 


2 
































melt migration in the mantle. In a review of existing attempts to model U-series melting 


models, Elliott and Spiegelman (2014) note that models generally reqnire a large set of 
assumptions such as melt-mineral partition coefficients, mantle porosity, source miner¬ 
alogy and petrology, the way in which melt is organized in the mantle, and the location 
and mechanism of chemical fractionation. Unsurprisingly, the combined outcome of 
these models is rather complex. Interpretations of ^^°Th disequilibria suggest that the 
mean flow speed between the melt source and ridge axis is around 1 m/yr, and shorter 
lived ^^®Ra disequilibria can be interpreted to require speeds >50 m/yr (Rubin and 


MacDougallj |1988t [Kelemen et al.| |1997t [McKenziel |2000t [Stracke et al.[ |2006t [Turner 


and Bourdon, 2010). These speeds are fast enough to require magma flow by some 


form of channelised melt transport. However, existing models of chanelized flow appear 


to be incapable of reproducing the ^^°Pb disequilibria signatures observed by Rubin 


et ah (2009), thus ^^°Pb dehcits are unexplained by current theories of channelized 


melt migration. 

How U-series disequilibria in MORE relate to chemical signatures of mantle hetero¬ 
geneity is the focus of only a small number of investigations. In a study of the Southeast 


Indian Ridge, Russo et ah (2009) investigated nonsystematic geographic variation of 
U-series disequilibria with axial depth, crustal thickness and ridge morphology. The 
authors explain the observations by along-axis variation in the melt supply, which they 
suggest is a response to long wavelength temperature gradients and small variations in 
source heterogeneity. Several other studies note a weak, positive correlation in noisy 
data between the amplitude of ^^‘^Th disequilibria and enrichment of the source in in¬ 


compatible elements (Lundstrom et ah, 1995 Lundstrom, 2000 Elliott and Spiegelman 


2003 Kokfelt et ah, 2003 Koorneef et ah, 2012). Since Th and U are fractionated more 


strongly in the presence of garnet, one explanation for this trend is that larger ampli¬ 
tude disequilibria result from a relatively larger proportion of melt originating in the 
garnet stability held (Lundstrom et ah, [1998 ). Yet, as Russo et ah (2009) demonstrate, 
it is non-trivial to make an explicit connection between the amplitude of disequilib- 
ria and magmas from garnet-bearing source rocks. Furthermore, it is unclear how the 
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dynamics and time scale of magma flow might contribnte to the observed signals. 

The cnrrent dynamical nnderstanding of melt migration is fonnded on mathematical 
theories and nnmerical models for the physics of magma flow in porons, compactible 


mantle rock (McKenzie, 1984; Fowler, 1985 Ribe, 1985 Bercovici et ah, 2001). A goal 
for many models of magma flow is to explore and nnderstand the dynamics of melt 
localization. Among these, channel formation is attribnted to the reaction infiltration 


instability (Ortoleva et al., 1987) that in the context of melt migration, describes a posi¬ 
tive feedback between magma flnx and dissolntion of porons mantle rock. This feedback 
causes a net transfer of mass from the solid mantle to the magma, locally increasing 


the porosity and melting rate, and results in the formation of channels (Kelemen, 1990 


Kelemen et ah, 1995 Aharonov et ah, 1995). As a consequence, models of this pro¬ 


cess must treat the physics and thermochemistry of magma flow as a coupled system. 
For simplicity, previous studies generally restrict the problem to magma flow in an 
upwelling column of partially melting mantle rock, neglect conservation of energy and 


parameterise melting using kinetic theory ( 

Spiegelman et al. 

2001 

Spiegelman and 

Kelemen, 2003 |Liang and Guo, 2003 Liang et ah, 2010; Schiemenz et ah, 2011 Hesse 

et al. 

2011 

)• 


Other models by Katz (2008, 2010), Hewitt (2010), Weatherley and Katz (2012) and 


Katz and Weatherley (2012) differ by adopting an energetically consistent approach to 


melting. These models capture the onset of partial melting and suggest that latent 
heat effects suppress the emergence of channels in the situation of magma flow through 


a chemically homogenous mantle. However, models by Weatherley and Katz (2012) 


and Katz and Weatherley (2012) show that channels emerge in a heterogeneous mantle, 
where spatial variations in composition are considered to be synonymous with variations 


in the fusibility of mantle rock. Weatherley and Katz (2012) shows that more fusible 
heterogeneities melt preferentially, and that the resulting flux of magma can be sufficient 


to nucleate high porosity channels. In a companion study, Katz and Weatherley (2012) 
considers the effect of different topologies of mantle heterogeneity and plate spreading on 
melt migration; they hnd that both heterogeneity and plate spreading have important 
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consequences for the arrangement of channels and delivery of melt to the ridge axis. 


In this contribution we extend the models of Katz and Weatherley (2012) to quan¬ 


tify the effect of source heterogeneity and spreading rate on the speed and time scale of 
melt migration. Section [^outlines the theory and numerical model of magma flow in 
heterogeneous mantle rock, summarises the petrological system (which is a key compo¬ 
nent of the model) and provides details on the approach used to determine the speed 
and time scale of melt migration. Solutions to the numerical experiments are presented 
in section Results are presented in dimensional form. Those presented at the start of 
the section aim to characterise the time scales and flow speeds predicted by the models, 
while those at the end of the section illustrate the effect of heterogeneity. The results 
are discussed in section where particular attention is paid to how the time scale of 
melt migration depends on spreading rate and mantle heterogeneity. A scaling analysis 
explores how the results might change under different parameter regimes; the assump¬ 
tions and limitations inherent in the models are also noted. The section concludes 
with a discussion on the implications of the results for our understanding of U-series 
disequilibria in MORE, including tentative prediction of their covariance with source 
composition. 

2. The model 


2.1. Model configuration 

The situation we consider is of magma flow, melting and melt-rock interaction 
in passively flowing, compositionally heterogeneous mantle beneath mid-ocean ridges 
(hgure [^. Subsolidus mantle rock flows into the bottom of the domain, upwells and 
melts adiabatically, and the melt then makes its way towards the surface where it is 


erupted or freezes into the lithosphere. Previous models by Katz (2008, 2010) and 


Katz and Weatherley (2012) form the basis for those presented here, although they are 


modihed in order to derive the predictions of the speed and time scale of melt migration 
beneath ridges. A brief summary of the behaviour predicted by the model is given in 
section o 
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Figure 1: A schematic illustration of the model setup. Lithospheric plates (light grey regions) spread 
apart at a known rate, driving large scale flow in the solid mantle (dark grey arrows). Decompression 
causes the mantle to melt within an unknown region indicated by the dashed line. Within this region 
solid rock and magma coexist and the melt segregates by porous flow (blue arrows) 


The models solve equations for conservation of mass, momentum, composition and 
energy in a thermally and chemically diffusive system with two phases and two thermo¬ 
dynamic components. In other words, they seek solutions for the coupled dynamics of 
solid mantle rock and liquid magma, and employ a basic petrological model to allow for 
melting and freezing. In the model framework mantle rock is considered to be a com- 
pactible, crystalline solid that deforms by creeping flow; magma is modelled as a low 
viscosity liquid that migrates by porous flow. In partially molten regions, the magma 


and mantle are assumed to coexist and interpenetrate (McKenzie, 1985 Fowler, 1985 


Ribe 

1985 

Bercovici et ah 

2001 

). The boundary conditions for the models are identi- 

cal to those used by 

Katz 1 

2010) 

, the only difference being that the internal boundary 


condition for melt extraction is pinned to the apex of the melting region, rather than to 
the boundary at the mid-line of the domain. Full details for this condition are provided 
in 


Weatherley (2013). An additional deviation between the models here and previous 


variants is that mantle flow is considered to be a passive response to plate spread¬ 
ing; variation in the gravitational body force on the two-phase aggregate is neglected. 
The equations are discretized using a hnite volume approach on a staggered mesh with 
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Figure 2: A binary solid solution with linearized solidus (lower surface, green lines) and liquidus (upper 
surface, black lines) constitutes the phase diagram. Co marks the bulk composition at zero pressure, 
Cf and Cm indicate the compositions of coexisting magma and rockm also at zero pressure. The mass 
fraction of the liquid is (Co — C/)/(Cm — Cf). 


1.25 km grid spacing in each direction. The governing eqnations are solved using the 


method described by Katz and Weatherley (2012). 


2.2. Phase relations and mantle heterogeneity 

Mass transfer between the solid and liquid by melting and freezing is an important 
feature of the model. To allow for this, we follow the examples of Ribe (|1985); Katz 


(2008, 2010); Hewitt (2010); Katz and Weatherley (2012) and Weatherley and Katz 


(2012) and assume that coexisting rock and magma are in local thermodynamic equi¬ 


librium. An important consequence of this assumption is that magma and mantle rock 
are in local chemical and thermal equilibrium in all partially molten regions. A simple 
phase diagram (hgure comprising two components of different fusibility describes the 
petrological system. At every time step the phase diagram is applied in each grid cell 
in order to relate the bulk composition, bulk enthalpy, and lithostatic pressure to the 
temperature, phase compositions, and rate of mass transfer between rock and magma. 
Since the phase diagram corresponds only to a hypothetical, highly simplihed petro- 
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Figure 3: The compositional heterogeneity fields used to initialize the models, shown as a perturbation 
to the solidus temperature, ATq = Ms{Cm — Cq). Cq — 0.5 for both fields. Dashed lines in (a) and 
(b) show the size of the computational domain relative to the heterogeneity fields. The heterogeneity 
fields are larger in the z-direction in order to evolve the bottom boundary condition on C as material is 
advected into the domain, (a) ‘Random blob’ style heterogeneity . 10% of the area is occupied by more 
fusible material, (b) ‘Smoothed noise’ type heterogeneity, (c) Directionally averaged power spectrum 
for the heterogeneity field in panel (b). It shows approximately uniform power for wavelengths > 11 
km and greatly reduced power for shorter wavelengths. 


logic system, compositional variations are expressed in terms of their perturbation to 
the solidus temperature relative to a reference composition, or equivalently referred to 
here as fusibility. This is expressed as ATq = M{Cm — Cq), where M is the slope 
dT/dC of the solidus. Cm is the local composition of the mantle rock and Cq is a 
reference composition. 

Figures and show the two contrasting heterogeneity helds used to initialize 
the models. They are constructed using stochastic algorithms that generate spatial 
variations in the ratio of the two thermodynamic components. Full details of the al¬ 
gorithms are provided in Katz and Weatherley (2012). In £gure|^, the ‘random blob’ 






































heterogeneity model consists of randomly shaped, isolated regions that are nniformly 
enriched in the more fnsible component. Prior to the onset of melting, the enriched 
regions pertnrb the local solidus temperature by —50 K, have a characteristic size of 
10 km and occupy an area fraction of 0.1. The composition of the non-enriched matrix 
dehnes Cq. The ‘smoothed noise’ model (£gure[^) is produced by allowing the ratio of 
the components to vary smoothly and continuously over all wavelengths greater than a 
specified cut-off. The spectrum has equal power at all wavelengths greater than 10 km 
and negligible power for wavelengths <5 km (£gure|^). Prior to the onset of melting, 
ATc varies by up to ±50 K and Cq is set to the 2-norm of the heterogeneity held. The 
inspiration for these conjectural descriptions of mantle heterogeneity lies in (a) recogni¬ 
tion that regions of the mantle enriched in ancient oceanic crust are more fusible than 


those that are not (Hirschmann and Stolper, 1996), and that the mantle contains re¬ 


fractory domains that are essentially unsampled by melting beneath oceanic volcanoes 


(Liu et ah, 2008), (b) uncertainty surrounding the lithology of mantle heterogeneity, 
and whether recycled oceanic crust melts and enriches the surrounding mantle prior to 


the main melting event beneath ridges (e.g. Yaxley and Green, 1998), (c) uncertainty in 


the topology of heterogeneity (Allegre and Turcotte, 1986 Kellogg and Turcotte, |1990 ), 
(d) analyses of scattered seismic waves that suggest chemically distinct heterogeneities 


with length scales of around 10 km exist in the mantle (Kaneshima and Helffrich, 1999 


2003), (e) geochemical estimates that the MORE source region contains around 10% 
recycled oceanic crust (e.g. Wood, 1979, Sobolev et al.| 2008; Shorttle and Maclennan 


2011), and (f) the need to represent mantle heterogeneity on scales larger than the grid 


resolution of the simulations (1.25 km). 


2.3. Determination of melt migration timescale 

To determine the speed and duration of melt migration predicted by the model, the 
domain is seeded with Lagrangian particles that track the motion of individual huid 
parcels. The magma velocity is given by 


K 

= ^m- — [VPrf ± Apg] , 


( 1 ) 
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where K = kocf)^ is permeability (fco is the permeability coefficient, 0 is porosity and n 
is the permeability exponent) /x is the viscosity of the magma, Pd is the dynamic flnid 
pressure, Ap is the density difference between the solid mantle and liquid magma and 
g is the acceleration due to gravity. Equation [T] states that the velocity of the magma 
depends on the matrix velocity, dynamic pressure gradients, and buoyancy forces that 
are modulated by permeability and the viscosity of the magma. The tracer particles 
are dehned to be perfectly incompatible; they are released from the solid at the onset 
of melting, can be re-incorporated into the solid by freezing, and can exit the domain 
at the ridge axis. The particles are accurately tracked by updating their positions at 
each time step in a two stage process. First the velocity held is linearly interpolated to 
the position of each particle. Then using the interpolated velocity, the position of each 
particle is updated over one half of a time step. The process is then repeated and the 
tracer particles are moved over the remaining half of the time step. 

The tracer particles are initialized by placing one in the centre of every grid cell. As 
the particles move through the domain, they record at every time step their velocity, 
position, the porosity and composition of the matrix within their host grid cell, and 
the model time. Those that exit the domain at the ridge axis are hltered to exclude 
any that {i) are not initialized in subsolidus mantle for which the degree of melting is 
zero, and {ii) are refrozen into the matrix at any point between the melt source and the 
ridge axis. Data from particles that pass through this hlter are used to compute the 
duration and speed of melt migration. In subsolidus regions, the motion of the tracer 
particles is equal to the solid mantle flow held. Since the how held approximates the 
corner how solution, the upwelling rate is fastest directly beneath the ridge axis and 


slower in the hanking regions (Batchelor, 1967). On long time scales a consequence 
of this is that the rate of delivery of tracer particles to the melting region is uneven 
across its width. To reduce any sampling bias the tracer particles are reinitialized 
periodically on a time scale of 350 kyr. Particles are not reinitialized more frequently, 
since magma migrating to the ridge on shorter time scales are capable of preserving 
230Th disequilibria generated at the base of the melting region. 
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3. Results 


3.1. Overview of dynamics 

This study is based around three suites of numerical experiments. One suite is 
initialized with a chemically homogenous mantle, another with the random blob het- 
eroegenity model (figure]^), and the final suite with the smoothed-noise heterogeneity 
model (figure [^). Each suite comprises four simulations with half-spreading rates Uq of 
2, 4, 6, and 8 cm/yr. Some key aspects of the model behaviour are outlined below, but 


the reader is referred to Katz and Weatherley (2012) for a complementary description 
and analysis. 

Figure]^ illustrates the basic dynamical behaviour of the system with a snapshot 
from a simulation run with no initial heterogeneity and a half spreading rate of 4 cm/yr. 
Models run with a homogenous mantle do not predict channelized flow; instead melt 
migrates by dynamically stable, diffuse porous flow. Melting is confined to a roughly 
triangular-shaped region extending about 75 km deep into the mantle and over 100 km 
to either side of the ridge axis. Throughout this region, coexisting magma and mantle 
rock are in local chemical and thermodynamic equilibrium (section |2.2[ ). Compaction 
maintains the porosity at low values, typically around < 1%. However, the porosity 
varies systematically throughout the melting region, largely as a response to the pattern 
of mantle upwelling; porosity is largest directly beneath the ridge axis where the rate of 
mantle upwelling (and decompression) is greatest, and decreases towards the flanks of 
the melting region. This porosity distribution sets the overal pattern of melt migration 
towards the ridge axis. Within c. 30 km either side of the spreading axis, melts 
are focused to the ridge axis by buoyancy and lateral pressure gradients, agreeing with 


predictions by Spiegelman and McKenzie (1987) and Yinting et ah (1991). At distances 


greater than 30 km, the mode of melt focusing resembles the model by Sparks and 


Parmentier (1991), in which magma rises vertically to the base of the lithosphere, and 


then uphill along the top of the melting region to the ridge axis. 

In numerical experiments initialised with a heterogeneous mantle, magma flow lo¬ 
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Figure 4: Snapshots from representative simulations run with contrasting styles of mantle heterogene¬ 
ity. In each case the half spreading rate is 4 cm/yr and the model time is 1.9 Ma. Panel (a) shows 
results from a simulation run with a homogenous mantle. Colours and solid dark grey contour lines 
show the porosity and extent of the melting region; the dotted streamlines show the solid mantle flow 
field; the solid blue-green lines plot the trajectories of some tracer particles that travel to the ridge 
axis. Panels (b)-(d) show results from a simulation with the random blob type heterogeneity, and 
panels (e)-(g) show results from a simulation run with the smoothed noise heterogeneity. Panels (b) 
and (e) express mantle heterogeneity as a perturbation to the solidus (grey line) temperatures, (c) and 
(f) show porosity, and d and g show temperature. The solid grey lines contour temperatures between 
0°C and 1300°C. The models are subject to the parameter regime outlined in table [ST| 
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calizes into high porosity channels. Fignre|^ shows snapshots of the matrix composition 
(fignre S>. e), porosity (fignre |^,f) and temperatnre ((hgnre ,g) from representative 
experiments with Uq = 4cm/yr. In these simnlations, the arrangement of channels 
and dynamics of focnsing strongly depend on the topology of mantle heterogeneity. In 
the random blob experiments (fignres S>. c, d) some channels are narrow and vertical, 
while others coincide with smeared-ont heterogeneities. Melt focnsing to the ridge axis 
is realised partly throngh channelized flow, and partly by flow along the lithosphere- 
asthenosphere bonndary. In the smoothed noise experiment (fignres |^, f, g) channels 
are wider and focns magma to the ridge axis more efficiently. Channels are deflected 
aronnd less fnsible regions of the mantle, and tend to coincide with compositions that 
he between the average mantle and those moderately enriched in the more fnsible com¬ 
ponent. Althongh coexisting magma and mantle rock are in local chemical and thermal 
eqnilibrinm in all partially molten regions, significant chemical diseqnilibrinm can exist 
between magmas within channels and magmas ontside of channels. This is becanse the 
timescale for melt transport throngh channels is mnch shorter than the timescale for 


diffusive equilibration of the channel and host mantle (Spiegelman and Kenyon, [1992 ). 

In these examples, variation of porosity within the melting region is a consequence of 
combined mechanical (compaction) and thermodynamic (melting and freezing) effects. 
Weatherley and Katz| (2012) showed that the melting rate T can be expressed as the 


sum of contributions from mantle upwelling, reactive flow and thermal diffusion. The 
melting rate is such that in a homogenous mantle, mantle upwelling accounts for over 
95% of the total melting rate. But in a heterogeneous mantle where flow is channelized, 
upwelling can account for a minor fraction of the total melting rate (as small as 20% 


in simulations by Weatherley and Katz (2012)), with the remained realised through 
reactive dissolution and thermal diffusion. Thus porosity variations in a heterogeneous 
mantle are largely a response to local thermodynamic effects. 
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3.2. Duration and distance of melt migration 

Figure plots the time taken for tracer particles to travel from the base of the 
melting region to the ridge axis as a function of the distance that they travel. The 
duration of melt migration for any tracer particle is also referred to as the travel time, 
r, and the distance travelled by tracer particles d is interchangeably referred to as the 
path length. These quantities are dehned as 

'T texit treleasej (2) 

d= |v/(t)|dt, (3) 

J release 

where treiease IS the time at which a particle is released for the hrst time from the solid 
by melting and texit is the time at which the same particle exits the melting region at 
the ridge axis. 

Figure shows results from experiments initialised with a homogeneous mantle; 
the travel time increases approximately linearly with the path length and is shorter 
when the spreading rate is fast. Path lengths range from 50-75 km when Uq = 2 cm/yr 
to 55-95 km when Uq = 8 cm/yr. For each experiment the range of travel times is 
approximately 150 kyr, and decreases from 200-350 kyr when Uo = 2 cm/yr to 90- 
240 kyr when Uq = 8 cm/yr. The arrays of data are more densely populated with 
tracer particles representing shorter path lengths and travel times since the upwelling 
rate is fastest and porosity is highest directly beneath the ridge axis. The spreading 
rate dependence of travel time is discussed more fully in section |4.1[ 

Figures It and show results from experiments run with the random blob and 
smoothed noise heterogeneity models. To leading order, these data show the same 
trends as experiments with a homogenous mantle (£g. |^); the travel time increases 
with path length and decreases with spreading rate. However, the data in hgures |^,c 
exhibit much scatter and the correlation between path length and travel time is weaker. 
Compared to results from the homogenous case, the data span a wider range of travel 
times and path lengths. For a half spreading rate of 2 cm/yr, the range of travel times 
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Figure 5: Travel time for individual particles plotted as a function of path length for half-spreading 
rates Uq 2, 4, 6, and 8 cm/yr. Panel (a) shows data from simulations with a homogenous mantle. 
Grey lines correspond to the magma flow speeds indicated. The number of particles reaching the ridge 
axis n in each simulation, arranged in order of increasing Uq is 2168, 6663, 11765 and 17821. Panel 
(b) shows data from simulations where the heterogeneity is manifested as random blobs. For these 
simulations n = 363, 805, 1380 and 1828. Panel (c) plots data from simulations run with smoothed 
noise heterogeneity. For these runs n = 30, 103, 278 and 575. 
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is approximately 50-350 kyr and path lengths span 25-80 km; for Uo = 8 cm/yr travel 
times span 25-300 kyr and path lengths range between 25 km and 110 km. The arrays of 
data in £gures|^,c are more densely populated for intermediate path lengths and travel 
times. They show that the duration of melt migration is shorter in the channelized flow 
regime by around 10-20%. The sensitivity of these results to the parameter regime for 
the model and discussion on their correspondence to U-series disequilbria is presented 
in sections 14.31 and 14.51 

3.3. Speed 

Figure shows normalized histograms of the particle mean speed for each of the 
numerical experiments presented in hgure The mean speed for particle p is dehned 
as (vp) = Tp/dp, where Tp and dp are the particle-specihc travel times and path lengths 
dehned in equations and Figures [^-d show data from experiments run with a 
homogenous mantle. These histograms have a strong negative skew, revealing that 
particles with the fastest speeds are the most abundant. From hgure it is evident 
that the fastest particles are generally associated with shorter travel times and path 
lengths. The mode of the binned data increases with spreading rate from 0.20-0.25 m/yr 
at t/o = 2 cm/yr to 0.65-0.70 m/yr when Uq = 8 cm/yr. The range of mean particle 
speeds also increases with spreading rate from 10 cm/yr when Uq = 2 cm/yr to 35 cm/yr 
when Uq = 8 cm/yr. 

Results in hgures[^-h and|^-l correspond to experiments run with the random blob 
and smoothed noise heterogeneity models. These histograms have variable, positive 
skew and show a factor of ~2 increase on the range of mean particle speeds predicted 
by experiments run without mantle heterogeneity. Mean particle speeds range from 
0-0.6 m/yr when Uq = 2 cm/yr to 0.15-1.3 m/yr when Uq = 8 cm/yr (hgures &-!)■ 
Generally, histograms for the smoothed noise cases have a more irregular structure than 
for those corresponding to experiments run with the random blob heterogeneity model. 
The mode of the binned, mean particle speeds is approximately independent of the 
presence or topology of mantle heterogeneity. 
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Figure 6: Normalized histograms showing the local mean speeds of particles. Each column shows results 
from simulations with the same half-spreading rate and each row presents data from simulations with 
different styles of heterogeneity as the initial condition. Bin widths are 5 cm/yr, n is the number of 
data points contained within the histogram, and s is is an estimate of the population skewness. 

3.4- Effect of composition 

Figure [^reveals the effect that the source fusibility at the on set of melting, ATc, 
has on the duration and speed of melt migration. Data are taken from experiments run 
with the smoothed noise heterogeneity model only. 

Figure shows that the path length correlates negatively with AT^. This neg¬ 
ative correlation is a consequence of the fact that more fusible rocks (more negative 
ATc) start to melt deeper. The combined effects of the shape of the melting region, 
channelized flow, and mantle heterogeneity mean that a range of path lengths can be 
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associated with any given composition. Data are distribnted relatively evenly across 
sonrce compositions between —20°C < ATc < 20°C bnt are less abnndant for more 
extreme valnes. Fignre shows that to leading order, the travel time and sonrce 
compositions are related by a weak to moderate negative correlation for half spreading 
rates of 4, 6, and 8 cm/yr, althongh no correlation is evident for Uq = 2 cm/yr. The 
leading order trend arises from more fnsible heterogeneities melting at deeper depth, 
while the weak correlation implies that the connectivity of channels and permeability 
strnctnre does not vary systematically with sonrce composition. 

Fignre plots the mean particle speed against ATq. Althongh no obvions trend 
relates the two qnantities, it can be seen that mean speeds are slower for regions where 
the fnsibility pertnrbations are large (|ATc| > 20°C) and faster where fnsibility is 
pertnrbed to a lesser extent. This resnlt reinforces the implication from fignre of 
non-systematic variations in permeability strnctnre and channel connectivity. 

To complement these data, the speed and composition for each particle were sam¬ 
pled at intervals of abont 1000 years after the onset of melting to give the instantaneons 
speed and instantaneons composition. These results, shown in figure [^, exhibit negli¬ 
gible correlation between speed and composition along trajectories of tracer particles. 
However, the greatest range and fastest speeds are observed in regions where mantle 
heterogeneity perturbs the solidus temperature between —15°C and -|-10°C. In these 
regions, instantaneous speeds are, in some instances, as fast as 30 m/yr, but are often 
as low as a few decimetres per year. For more extreme compositions, the range and 
magnitude of instantaneous speeds is much reduced. In these regions the maximum 
magma speeds are around 2 m/yr, but the vast majority of tracers passing through 
more highly enriched and depleted regions have speeds of a few centimetres to a few 
tens of centimetres per year. 
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Figure 7: Effect of composition on the time scale and speed of melt migration for experiments run 
with the smoothed noise heterogeneity model, (a) Variation of travel time with source composition 
AT. Black lines show least squares best fit lines for the different simulations and are labeled as Uj;, 
where x is the half spreading rate in cm/yr. The corresponding sample correlation r for each of the 
best fit lines is given in the upper right of the figure, and Ux is the number of data points from each 
simulation, (b) Path length as a function of source composition, with black lines as for panel (a). Panel 
(c) plots mean particle speed as a function of source composition; panel (d) shows variation of the 
instantaneous speed with instantaneous composition, ATc^nst.) expressed in terms of a perturbation 
in solidus temperature. The sampling interval is approximately 1000 years. 


19 





















4. Discussion 


4.I. Effects of plate spreading rate 

Two leading-order trends emerge in the results: that shorter melt-transport path 
lengths and faster spreading rates are both associated with shorter durations of melt 
migration. These trends arise because plate spreading conhnes melting to a roughly 
triangular-shaped region of non-constant porosity. Analysis of the governing equations 
in ID by 


Hewitt 


(2010) reveals that the magma speed w scales with where 


W is the upwelling rate of the matrix and n is the permeability exponent (Bear, 1972 


Turcotte and Schubert, 2002 Costa, 2006). For mantle convection patterns resembling 


corner flow, the half spreading rate dehnes a reference upwelling rate Wo = 3/2f/o 
(Batchelor, 1967), thus w oc Figure i shows that this scaling accounts for 

the spreading rate dependance of the travel time in a homogenous mantle. In the 
experiment with Uq = 2 cm/yr, tracers associated with short path lengths {d < 65 km) 
have anomalously slow speeds. This is a consequence of the slow upwelling causing 
smaller melting rates, lower porosity, lower permeability, and hence slower segregation. 


4 . 2 . Effects of heterogeneity 

Compared to plate spreading, mantle heterogeneity has a more complicated effect 
on melt migration. While the fact that less fusible mantle melts at shallower depths 
readily explains the negative correlation between travel time and source composition, 
the relationship between speed and composition is more difficult to understand. To help 
explore why faster speeds are associated with source and regions where perturbations to 
the fusibility are relatively minor (|ATc| < 20°C, hguresj^, d), hgurecompares the 
spatial covariation of mantle heterogeneity with porosity for a representative experiment 
run with Uq = 6 cm/yr and the smoothed-noise heterogeneity model. In panels a and 
c, black lines contour porosities > 1% to roughly show the location of channels and 
magma ponds. In this experiment, channels are conhned, more or less, to the region 
within 100 km of the spreading axis and depths < 60 km. The channels accommodate 
the fastest magma speeds (hgure [^) and broadly coincide with mantle of intermediate 


20 


















3 



Normalized path length 


Figure 8: Scaled travel time versus nondimensional distance for the experiments conducted with no 
mantle heterogeneity. Travel time r is scaled by t* = (3/2{7o)^“^/"'/Z, where I is a characteristic length 
scale. 

fusibility (—15 K < ATc < +10 K); in regions of more extreme fusibility the flow 
speeds are slower. This correlation is also evident in hgure [^. 

This relationship between composition and magma speed is a consequence of the 
energetics of melting. The imposition of local thermodynamic equilibrium pins the 
temperature of partially molten regions to the solidus temperature. Thus heat diffuses 
away from locally depleted regions into the surrounding more fusible and cooler mantle. 
In the most depleted regions {ATc > lOK), therefore, melting is suppressed, porosity is 
destroyed by compaction and the resulting small permeability accommodates slow in¬ 
stantaneous magma speeds and acts to decrease the mean speeds of magmas originating 
in these regions. Flow speeds in the most enriched regions {ATc A — 15K) are slow for 
a different reason: enriched heterogeneities start to melt before the surrounding man¬ 
tle and are surrounded by a halo of impermeable mantle rock. This prevents magma 
from flowing into the channel network; within the heterogeneity the magma speed |vj| 
approximates the mantle speed jv^l, which results in slow instantaneous and mean 
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Figure 9: Covariation of porosity and composition from a representative experiment initialized with 
the smoothed noise heterogeneity model and Uq = 6 cm/yr. (a) Colours show composition in terms 
of a perturbation to the solidus temperature; grey lines contour (p = 1%. (b) Colour shows porosity; 
grey solid lines contour ATq = —15°C; grey dotted lines contour ATq = +10°C. (c) Covariation of 
large porosities and regions of enriched and depleted mantle. Black lines contour (j) = 1%, blue lines 
contour ATq < —15°C in 5°C intervals, and red^^nes contour ATc > +10°C in 5°C intervals, (d) 
Magma speed |v/| with grey lines showing the 0.5% porosity contour. 



















speeds. It is interesting to note that, since the average flow speed for each particle is, in 
part, a response to local thermodynamic conditions, different heterogeneity configura¬ 
tions could result in speed-composition distributions different to those shown in figures 
[H:: and d. 


4-3. Effects of model parameters 

An important point about the results above is that they are sensitive to the model 


parameters, and particularly those that directly affect the magma velocity. Appendix 


[^presents a nondimensionalization of equationfor vj and shows that the nondimen- 
sionalizing factor is 7^ = kQApg/fi, where is the reference permeability, Ap is the 
density difference between the matrix and magma, g is the acceleration due to gravity 
and p is the viscosity of the magma. Of these, g is accurately known, possible values 


for Ap vary to within a few percent of the value given 

surements of p span one order of magnitude (1-10 Pa-s, 

in table 

Dingwell 

Al (500 kg m^), mea- 

(1995)) and estimates 

of ko lie in the range 10 ® m^ to 10 ^ m^ ( 

Paul 

,11997; 

Wark and Watson 

1998 

Faul 


2001 Wark, 2003 Connolly et ah, 2009). Uncertainty in these parameters means that 


the range of plausible values for TZ spans some 6 orders of magnitude. For comparison, 
a permeability constant of /cq = 10“^m^ was selected for the model runs above, and the 
corresponding value of TZ lies approximately in the centre of the possible range. 

To investigate how the mean particle speed (vp) and travel time r vary with TZ, 
the model was reconfigured to consider the dynamically simpler case of a column of 
upwelling, chemically homogenous mantle rock. Figure [l0| plots the predicted (vp) and 
r for values of TZ that span 5 orders of magnitude. The data are normalized by results 
from the model run with the parameter regime specified in table |A1[ The figure shows 


that (vp) scales with TZ^^'^ (figure 10 r), where n is the permeability exponent and that 
T varies with This scaling agrees with analysis of the governing equations in 


one dimension by Hewitt (2010); it indicates that an change in 7?. by a factor of 10'* 
will yield an order of magnitude change in (vp) and r. Provided that the scaling holds 
for models with a heterogeneous mantle, setting TZ = 10^7^* would reduce the average 
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Table Al: List of parameters and preferred values. 


Symbol 

Meaning 

Value 

K 

Thermal diffusivity 

1 X 10“® m/s^ 

L 

Latent heat of fusion 

5.5 X lO^J/kg 

Ms 

Slope dT jdC of the solidus 

200 K 

Ml 

Slope dT/dC of the liquidus 

200 K 

T 

Potential temperature 

1696 K 

To 

Solidus temperature for reference 

mantle at 0 kbar 

1605 K 

Co 

Reference mantle composition 

0.5 

AC 

Difference in composition between 

rock and magma in thermodynamic 

equilibrium 

0.1 

AT 

Temperature scale 

MAC 

Ap 

Matrix-melt density difference 

500 kg/m^ 

a 

Coefficient of thermal expansion 

3 X 10“^ K-i 

Cp 

Specific heat capacity 

1200 J/kg/K 

7“^ 

Clapeyron slope 

60 K/GPa 

V 

Chemical diffusivity 

1 X 10“® m/s^ 

9 

Acceleration due to gravity 

9.8 m/s^ 

ko 

Reference permeability 

1 X 10“^ m^ 

7max 

Cut-off shear viscosity, maximum 

lO^spa-s 

Ax 

Grid spacing in x-direction 

1.25 km 

Az 

Grid spacing in ; 2 -direction 

1.25 km 


travel times to order 10"^ years and increase mean flow speeds to order 10 m/yr. 
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Figure 10: The effect of 7^ = koApg/^ on melt migration computed using the upwelling column model 
outlined in Weatherley and Katz| (2012|. The upwelling rate is 6 cm/yr and the grid spacing is 1 km. 
Other parameter values are given in table A1 (a) Effect of TZ on the mean particle speed (vp). TZ* and 


(vp)* indicate normalizing values computed from an experiment conducted with the parameter regime 


in table [Al] (b) Effect of TZ on the duration of melt migration r, where r* is associated with the TZ* 
parameter regime. 


4.4- Limitations and assumptions 

The results above are also subject to a set of limitations arising from the assumptions 
that underpin the models and theory. An important assumption is that melt migration 


is accommodated entirely by porous flow (McKenzie, 1984 Fowler, 1985; Ribe, 1985 


Bercovici et ah, 2001[ ). Whilst a mass of geochemical and geophysical observations 


require the mantle to be porous and melt to migrate by porous flow (e.g. Johnson 


et al. 

1990 

Johnson and Dick 

1992 

The MELT Seismic Team, 1998; |Sims et al. 

1999 

), several studies also suggest that melt can flow through cracks ( 

Sleep, 

1988 


Fowler and Scott, 1996 Ito and Martel, 2002 Maaloe, 2003). This mode of magma 


flow is less well studied than reactive channelisation, thus uncertainty surrounds how 
widely it applies to the sub-ridge melting region. Nevertheless, magma flow through 
open cracks can accommodate much faster flow speeds than melt migration through 
high porosity channels; hence this mode of localisation lends itself well to arguments 


for preservation of deep-origin ^^°Pb disequilibria in zero-age MORE (Rubin et al. 


2005). Future models should seek to include a brittle fracture criterion in order to 
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assess the implications of fracture for melt extraction. 

An additional limitation lies in the two-dimensional description of the dynamics at 
mid-ocean ridges. This imposes an inhnite, out-of-plane length on all heterogeneities 
and melt migration features in the model and restricts mass and heat transport to be in 
the plane of the solution. Katz and Weatherley| (2012) note that this latter restriction 
could reduce the efficiency of channel formation. For instance, if a channel has a cylin¬ 
drical geometry in three dimensions, mass and heat can be supplied from around its 
circumference. Radial heat and mass transfer could potentially reinforce the channel¬ 
forming instability, resulting in higher porosities and faster magma flow speeds. Three 
dimensional consideration of melt migration is also important in the vicinity of ridge 
offsets. Previous studies show that hot mantle beneath ridge segment ends is cooled 


by the adjacent, older and cooler lithosphere (Magde and Sparks, 1997 Weatherley 


and Katz, 2010 Herbert and Montesi, 2011). This cooling effect is likely to suppress 


porosity generation, result in lower magma flow speeds, and focus melt towards segment 
centres. Exactly how the presence of mantle heterogeneity might modify melting and 
melt transport close to segment ends is not yet understood, but could have particu¬ 
lar relevance to the interpretation of bull’s eye gravity anomaly patterns along highly 


I'l.s.s 

Lin et al.| 

1990; 

Detrick 


et ah, 1995; Magde and Sparks, 1997). 


The petrological model is the most signihcant simplihcation in our theoretical frame¬ 
work for magma flow. Motivated by experimental studies of mantle rocks, the model 
captures variation of the mantle solidus with composition and pressure, and takes an 
energetically consistent approach to melting by decompression and reaction. These 
properties approximately align the model with natural system. For simplicity, however, 
the model system comprises two phases and two thermodynamic components. Whilst 
this precludes consideration of additional important petrological features of the sub¬ 
ridge mantle, such as variations in modal mineralogy, volatiles, phase specihc enthalpies 
of fusion, and reactive precipitation of pyroxene, it lends itself well to elucidating the 
dynamical consequences of mantle heterogeneity for melt migration. A particularly im- 
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portant result of the counterpart models by Weatherley and Katz (2012) and Katz and 


Weatherley (2012) is that the arrangement and existence of melt migration pathways 


is a response to local thermodynamic effects. This hnding is expected to hold both on 
smaller length scales and in simulations with an expanded petrological model; thus the 
current model provides a reasonable, first order description of the large scale dynami¬ 
cal behavionr. However, a more detailed consideration of the petrology may offer new 
insights to trace element redistribntion and generation of U-series diseqnilibria in the 
melting region. 


Similar to earlier models by Katz (2010), magma is extracted at the ridge axis. 
The internal bonndary that facilitates melt extraction is carefnlly tuned so that on 
time scales that are long compared to a time step, the rate of magma production 
approximately matches the combined rates of melt extraction at the ridge, freezing at 
the base of the lithosphere and exit of magma throngh the side walls of the domain. 
One consequence of this internal bonndary is a slight build-np of magma in the top 20 
km of the mantle beneath the ridge axis. In this region, the dynamic pressure gradients 
are such that the magma velocity is small, despite the large porosity. Regions of high 
porosity are also predicted to occur in isolated patches in the flanks of the melting 
region. Since the permeability increases as a power-law of porosity, magma within 
these regions can convect at speeds of hnndreds of metres per year, drastically redncing 
the size of the time-step in the simulations. To avoid this, the permeability is capped 
at a value corresponding to a porosity of 10%. For the most part, the permeability 
cap has negligible effect on the speed of melt migration within channels. It also has 
negligible effect on magma speeds in the region of slightly higher porosity immediately 
beneath the ridge axis, owing to the setup of the internal bonndary condition for melt 
extraction. 

An additional limitation of the model is the grid resolution. In the simulations above, 
the grid resolntion is 1.25 km, which limits the smallest length scale of heterogeneity 
and width of melt localization featnres. A key resnlt of the connterpart models by 


Weatherley and Katz (2012) and Katz and Weatherley (2012) is that the arrangement 
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and existence of melt migration pathways is response to local thermodynamic effects. 
Therefore, it is reasonable to expect a similar basic model behavionr in simulations 
run with much hner grid resolution. Alongside offering opportunities to explore the 
coupling between heterogeneity on sub-kilometer scales and magma flow speeds, hner 
resolution models offer greater possibilities for investigating how the number density of 
fusible and refractory heterogeneities inhuence the large scale pattern of melt migra¬ 
tion. Many studies of ophiolites and peridotite massifs demonstrate that mineralogical 


heterogeneity can be present on scales upwards of a few centimetres (e.g. Boudier and 


Coleman[|1981t [Gregory! [T984l|Reisberg et al.]|1991t|Bhchert-Toft et al.||l999t|Le Roux 


et al., 2007|). At these small scales, physical attributes of the mantle such as grain 


size 


variations and minor shear zones are likely to have consequences for melt localization, 
in addition to those from composition heterogeneity. Our understanding of how these 
different variables inhuence the dynamics of magma how and hence the time scale of 
melt migration will beneht from future studies that combine lithological and structural 
mapping of mantle peridotites on a range of scales with higher resolution numerical 
models with more advanced rheological laws than that used here. 


4-5. Geochemical implications 

A major scientihc challenge at present is to understand the chemistry of oceanic 
basalts in terms of source composition and processes that ahect magmas between the 
point of origin and eruption. A key result from the models above is that under conditions 
of strong mantle heterogeneity, melts migrate to the ridge axis on a broad range of time 
scales. We postulate that (i) this distribution in travel times maps on to signihcant 
natural variability in the U-series disequilibria observed along MORs, and (ii) that 
the results offer an explanation for the decoupling that is often observed between U- 
series disequilibria and other geochemical data, even in settings where there is evidence 


for mineralogical source heterogeneity (e.g. 

Kokfelt et al., 

2003 

Bourdon et al. 

2005 

Stracke et al. 

2003 

Prytulak and Elliott, 2009 

Russo et al. 

2009 

Koorneef et al. 

2012) 


To further understand the implications of these results, future studies should aim to 















































































constrain how properties of mantle heterogeneity and the melt transport network are 
linked. Possible avennes to investigate inclnde how the nnmber density, fnsibility and 
connectedness of heterogeneity relate to the topology and tortuosity of melt channels, 
whether channel tortuosity on smaller scales than that examined here has a signihcant 
effect on melt transport time, and how heterogeneity and magma flow together affect 
the rate of magma supply to ridges, and by extension their morphology. 

Within the context of the model results, it is also important to consider other 
processes that can have notable implications for the geochemistry of oceanic basalts. 
Among the most signihcant is mixing. Analyses of whole rock and melt inclusion 
chemistry by Maclennan et ah] (2003) and Maclennan (2008) demonstrate that erupted 
magmas are mixtures of many small batches of melt, with up to 30 batches of melt 


represented in a single hand specimen. Additionally, Rudge et ah (2013) showed that 
the mixing process is nonuniform with depth, with the deepest melts appearing to be 
more homogenized than shallower melts. Interestingly, they also found that compo¬ 
sitional arrays produced during mixing do not necessarily point towards the isotopic 
compositions of the source materials. A second set of processes important to the the 
preservation or homogenisation of primary melt compositions are those that occur in 


the crust and uppermost mantle. Rubin and Sinton (2007) showed that high melt sup¬ 
ply rates (typically associated with fast spreading ridges) promote larger and shallower 
accumulations of cooler magma, longer residence times in sub-ridge magma reservoirs, 
and consequently greater extents of homogenization and differentiation. In contrast, 
lower melt supply rates promote episodic magma accumulation in hotter, more poorly 
connected magma reservoirs, and erupted products that show high degrees of chemi¬ 
cal variability and lower degrees of differentiation. In the shallow crust, crystallisation 
processes and melt-rock reaction can obfuscate the relationship between primary melt 


chemistry and that of erupted magmas (e.g. Van Orman et ah , 2006; Maclennan, 2008). 


However, Rubin et ah (2009) note that few sections of the MOR system are sampled 
at high enough resolution to quantify how these processes, and by extension, results of 
numerical models, relate to MORE compositions and MOR structure. 
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Despite these additional complexities and limitations, the results above also have 
several general implications for our understanding of how theories of magma flow re¬ 
late to observations of U-series disequilibria in MORE. Under the parameter regime 


outlined in table the results indicate that magma flow in a heterogeneous mantle 
is most commonly fast enough to preserve ^^°Th excesses generated deep in the melt¬ 
ing region. Some particles travel slowly enough for ^^°Th initial excesses to decay to 
zero. However, signihcant excesses for these particles could be generated by ingrowth, 
or alternatively, their lack of ^^°Th disequilibrium could be obscured by mixing. For 
the same simulations, magma flow too slow to preserve deep-origin ^^®Ra disequilibria, 
yet section |4.3| suggests that under less moderate parameter regimes, the models will 
predict flow speeds sufficiently fast to enable attribution of some ^^®Ra disequilibria to 
a deep source. Under even the most favourable parameter regime, the current theoret¬ 
ical framework for porous melt migration is incapable of preserving ^^°Pb disequilibria 
generated within the melting region. 

To clarify this situation, and to address whether current theory can account for ^^®Ra 
disequilibria, future studies must consider the full range of disequilibrium-generating 
processes. These include chromatographic fractionation, ingrowth, degassing and diffu¬ 
sive partitioning between melt and cumulate rocks. Furthermore, the models must also 
be able to account for covariation of different parent-daughter nuclide pairs. Models 


by dull et ah (2002) suggest that magmas from different portions of the melting region 


must mix to produce the observed signature, and Elliott and Spiegelman (2003) demon¬ 
strate that so called ‘full reactive transport models’, which couple equations for element 
concentration to dynamical models similar to those described above, also predict an in¬ 
verse relationship between ^^®Ra and ^^°Th disequilibria, but fail to match the observed 
amplitudes. However, future modelling studies that consider U-series disequilibria in 
the context of mantle heterogeneity and the physics of melt migration may yield a more 
conclusive view on the presence of short-lived U-series disequilibria in oceanic basalts. 


30 










5. Conclusions 


Results from energetically consistent models of magma flow show that plate spread¬ 
ing and mantle heterogeneity have important consequences for the time scale of melt 
migration beneath ridges. Heterogeneity perturbs the thermal structure of the mantle 
that is established by plate spreading, and causes magma flow to localize into high 
porosity channels. Magma flows relatively rapidly through channels directly beneath 
the ridge axis where the upwelling rate of the porous mantle is greatest, and more slowly 
through channels on the flanks of the melting region. Flow speeds also increase with 
spreading rate, consistent with simpler models. The main effect of mantle heterogene¬ 
ity on the time scale of magma flow is realised through the formation of high porosity 
channels. However, magma flow is generally slower in more enriched and more depleted 
regions, and faster in regions of close to average fusibility. 

Our hndings suggest that the mean speed of channelized flow in a heterogeneous 
mantle is approximately 20% faster than diffuse flow in a homogenous mantle. For 
the parameter regime and spreading rates considered, individual packets of magma in 
a single simulation migrate to the ridge axis on time scales as short as 25 kyr and 
as long as 350 kyr, with mean speeds between 0.2 and 1.3 m/yr. Scaling analysis 
suggests that for plausible parameter regimes that capture a more permeable mantle, 
the duration of melt migration might be reduced by a factor of 10. These speeds 
are sufficient for mineral-melt partitioning deep in the melting region to contribute to 
^^°Th and, in some circumstances even ^^®Ra disequilibria observed in MORE. However, 
the models’ inability to explain observations of ^^“^Pb dehcits in MORE suggests that 
current theories of porous melt migration do not comprehensively describe the physics 
of magma flow. 

An important result is that mantle heterogeneity induces signihcant natural vari¬ 
ability in U-series disequilibria preserved in MORE, and thus the models offer one 
explanation as to why observations are often decoupled from other geochemical sig¬ 
nals and properties of ridge axes, even in situations where there is good evidence for 
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mineralogical heterogeneity in the mantle. 
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Appendix A. Parameter Values 


The parameter values used for the numerical experiments above are listed in table 
All below. 


Appendix B. Nondimensionalization 


To nondimensionalize equation the following scales are employed: 


(v/,Vm) = W^0(v/,Vn,), 


(B.l) 


K' = koK, 

(B.2) 

v' = —V 

H ’ 

(B.3) 

Pd = ApgHPd, 

(B.4) 


where nondimensional variables are indicated with a prime, H is the height of the 
domain and ko is the reference permeability. Application of these scales reveals that 
equation is nondimensionalized by the factor TZ = koApg/^. 
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